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Magnetohydrodynamic (MHD) turbulence in the solar wind is observed to show the spectral 
behavior of classical Kolmogorov fluid turbulence over an inertial subrange and departures from this 
at short wavelengths, where energy should be dissipated. Here we present the first measurements 
of the electric field fluctuation spectrum over the inertial and dissipative wavenumber ranges in a 
(3 > 1 plasma. The k~ 5 ^ 3 inertial subrange is observed and agrees strikingly with the magnetic 
fluctuation spectrum; the wave phase speed in this regime is shown to be consistent with the Alfven 
speed. At smaller wavelengths kpi > 1 the electric spectrum is softer and is consistent with the 
expected dispersion relation of short- wavelength kinetic Alfven waves. Kinetic Alfven waves damp 
on the solar wind ions and electrons and may act to isotropize them. This effect may explain the 
fluid-like nature of the solar wind. 

PACS numbers: 52.35.Ra Plasma turbulence; 52.35.Bj Magnetohydrodynamic waves 



Turbulence is ubiquitous in astrophysical plasmas; tur- 
bulent processes are thought to play a role in cosmic 
ray and energetic particle acceleration and scattering 
0, advection dominated accretion flows, and perhaps 
solar/stellar wind acceleration. Yet key aspects of the 
physics of turbulence in magnetized plasmas are poorly 
understood, in particular the physics of dissipation at 
small scales. The classical scenario of magnetohydrody- 
namic turbulence is thus: fluctuations in the plasma are 
driven at some large 'outer' scale and decay by interact- 
ing locally in k— space. Eddies at some scale k exchange 
energy with eddies at nearby spatial scales, possibly as a 
three- wave or higher order interaction [ToL [ill ll2T ] , with 
the resulting net flow of energy to smaller spatial scales 
(larger fc); this cascade of energy occurs over an 'iner- 
tial subrange' of k— space and can be shown to predict a 
power spectrum that scales as fc~ 5 / 3 . At the smaller scale 
of the ion thermal gyroradius, kpi > 1, the ions become 
demagnetized and the plasma can no longer behave as a 
simple fluid; the turbulent energy is then thought to be 
damped on the thermal plasma by Landau or transit-time 
damping. However, the details of this damping process 
are not known and there a few reported measurements 
in this regime of k— space. 

Observations of the magnetic spectrum show break- 
points at near kpi rs 1, above which the spectrum typ- 
ically becomes harder [^J [lj, [lj| . This has been inter- 
preted variously as evidence of kinetic Alfven waves , 
whistler wave dispersion|l7|. and ion cyclotron damping 
of Alfven waves [l|| ■ 

Here we report the first measured power spectrum of 
electric fluctuations in solar wind turbulence. The iner- 
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tial subrange is clearly evident and follows the magnetic 
fluctuation spectrum. At large wavenumbers kpi > 1, 
the electric spectrum is enhanced. 

Data are used from experiments on the Cluster space- 
craft. Cluster flies four spacecraft, as a controlled tetra- 
hedron, in an inclined orbit with apogee at 19 Earth radii 
(Re)- From December to May each year, the spacecraft 
exit the terrestrial magnetosphere on the dayside and 
make measurements in the solar wind. We use approxi- 
mately 195 minutes of data during the interval 00:07:00 
- 03:21:51 on February 19, 2002, when Cluster was at 
apogee and spent several hours in the ambient solar wind; 
all of our data is from Cluster spacecraft 4. 

The electric field is measured by the Electric Field 
and Waves experiment (EFW) experiment EFW is 
a double-probe experiment which measures the floating 
voltage of 8 cm diameter current-biased spheres extended 
on 44m wires in quadrature. These spheres, as well as 
the spacecraft, are illuminated by the Sun and emit pho- 
toelectrons which cause the surfaces to charge positive 
with respect to the plasma. The surfaces attract a return 
current of thermal electrons which provide the electrical 
coupling to the plasma. Systematic variations in this 
coupling, due to changing illumination or variations in 
surface properties and work function, are a large source 
of background noise in EFW at the spacecraft spin-period 
(4 seconds) and harmonics. This is discussed more below. 
EFW measures the electric field on two orthogonal sen- 
sor pairs in the spacecraft spin-plane at 25 samples/sec. 
These two components are rotated into X and Y com- 
ponents in the GSE (geocentric solar ecliptic) coordinate 
system. Since the GSE Y direction represents the ori- 
entation with best symmetry for solar illumination, this 
component of the electric field is generally less noisy; we 
use the GSE Y electric field E y for all of our analysis. 
However, at any given instance, E y is composed of data 
from all four electric field probes, each with slightly dif- 
ferent photocoupling to the plasma. We therefore apply 
a finite impulse response (FIR) filter to the data to notch 
out the primary perturbations at the spin-tone and some 
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harmonics. 

The magnetic field is measured by the FGM instrument 
; three-component magnetic field vectors are sampled 
at 22 samples/sec (SPS). In our analysis, we use the 
GSE Z component of the magnetic field B z for reasons 
that are explained below. Moments of the solar wind 
ion distribution (velocity density, and temperature) are 
computed from the ion spectrum measured by the CIS 
experiment^) . 

Figure 1 shows an overview of the data used in the 
following analysis; panels a) and b) are wavelet spectro- 
grams and will be discussed below. Panel c) shows the 
two components of measured electric field E x and E y 
in GSE coordinates. Panel (d) show the magnetic field 
data. Panels e), f), and g) show the plasma ion density, 
plasma ion /3$ (ratio of plasma to magnetic pressure), and 
Alfven Mach number. The average ion beta is f3i ~ 5, 
average Alfven speed va ~ 40 km/s, and the average 
solar wind velocity is v sw « (-347, 4.9, -32.6) km/s (in 
GSE coordinates), over the entire interval. During the 
interval between 00:30 and 00:50, the magnetic field is 
nearly tangent to the Earth's bow shock (as per a calcu- 
lation assuming straight field lines @); however, Cluster 
summary plots of electron and plasma wave data show 
no evidence of connection to the shock. All of our data 
is ambient solar wind. 

To compute power spectra, the electric field data E y 
(25 samples/sec) were subsampled onto the time tags of 
the magnetic field data B z (22 samples/sec) by linear in- 
terpolation; a total of exactly 2 18 points are used. The 
power spectral density (PSD) was computed using both 
Fast Fourier Transform (FFT) and Morelet wavelet 
schemes. The FFT was computed as follows: the data in- 
terval was divided into 64 contiguous ensembles of length 
4096 (182 seconds); this gives an inherent bandwidth of 
A/ w 1/186 Hz. To minimize spectral leakage, each 
ensemble was 'whitened' by applying a first-order differ- 
ence algorithm, the PSD was computed by FFT, then the 
spectrum was postdarkened |l| and divided by the band- 
width of the FFT. Since the data is prewhitened, no win- 
dow function was applied before the FFT. The electric 
field spectra were then 'cleaned' by interpolating over the 
narrowband spikes resulting from the spin-associated sig- 
nals described above. A final spectrum was computed as 
the average of the 64 ensembles. Figure 2 shows the FFT 
power spectra of E y and B z (in black) . Wavelet spectra 
were computed by first producing the (complex) FFT of 
E y and applying the spectral cleaning (interpolation) to 
the real and imaginary parts, at positive and negative 
frequencies. An inverse FFT restores the 'cleaned' signal 
and a Morelet wavelet spectrogram was computed from 
this cleaned E y , as well as the original B z . The wavelet 
has 136 log-spaced frequencies; the final wavelet PSD is 
computed as the square of the spectrum averaged over 
time. The wavelet PSD is also shown in Figure 2 (in 
red). The wavelet spectrum extends to lower frequencies 
than the FFT, which is composed of ensembles of smaller 
data intervals; however, these very low frequencies lie be- 



low the 'cone of influence' and are unreliable [9j. Here we 
restrict our interpretation to the region where the FFT 
and wavelet spectra agree. The FFT electric spectrum in 
Figure 2 shows clearly the effect of the notch filters and 
residual spin-harmonic spikes. The wavelet PSD, with its 
much larger bandwidth, mostly averages over these resid- 
ual features although a depression near the notched por- 
tion of the spectrum can be seen. The FFT and wavelet 
PSD spectra agree remarkably well for both electric and 
magnetic fields. 

Of course, our (human) scheme of measuring time 
means little to the solar wind plasma, so there is little 
reason to expect the data to be inherently organized by a 
power spectrum in Hertz. Since the solar wind is super- 
Alfvcnic (Figure 1), the phase speed va of the Alfvenic 
fluctuations is much less than the wind speed itself; hence 
the measured frequency spectrum is actually a Doppler- 
shifted wavenumber spectrum lu w kv sw . This is often 
called Taylor's hypothesis and might not be considered 
to hold at large wavenumbers, especially if waves are 
present with phase speeds greater than the solar wind 
speed (such as whistler waves). 

As discussed above, it is considered that the fluid-like 
behavior of the wind breaks down at near kpi w 1, there- 
fore kpi is a natural parameter for organization of the 
power spectrum. The top panel of Figure 3 shows the 
FFT and wavelet power spectra organized by kpi, instead 
of frequency. For the FFT spectrum, the local values of 
Ti, and \B\ are used to compute k = u>/v sw and 
the thermal ion gyror&dius p i — v^j 

averaged over 

each (186 sec) ensemble; the E y and B z power spectra 
are then interpolated onto a linearly-spaced set of val- 
ues kpi G (0.006, 10). Since solar wind parameters vary 
slightly in each ensemble, this also has the effect of smear- 
ing (averaging) over the narrowband interference in the 
FFT PSD of E y . The wavelet spectrograms are time- 
averaged on to 4 second intervals and then interpolated 
onto a set of log-spaced values of kpi\ panels a) and b) 
of Figure 1 show these scaled spectrograms as a function 
of time. In Figure 1 a) and b), the fluctuation power 
has been divided by fc~ 5 / 3 to highlight fluctuations above 
the average spectrum of the inertial range. The electric 
and magnetic wavelet spectrograms are then averaged to 
compute the composite spectra in panel a) of Figure 3. 

Between kpi 0.015 and 0.45, the wavelet and FFT 
spectra of electric and magnetic fluctuations show power 
law behavior with indices of fc -1 ' 7 , which is consistent 
with the Kolmogorov value of 5/3. Both SE V and 5B Z 
show breakpoints at near kpi « 0.45; the magnetic spec- 
trum becomes harder with a index k~ 212 , while the elec- 
tric spectrum becomes softer. As discussed above, hard 
magnetic spectra have been observed previously |14L Il5| . 
Above kpi « 0.45, the electric spectrum is power law like 
k~ 126 to kpi « 2.5. Above this second breakpoint, a ex- 
ponential exp (— fcpi/12.5) better fits the spectrum. At 
these higher wavenumbers, the electric field data is noisy 
and shows harmonics of the spin tone (as shown above) . 
To test the validity of this data, we perform two analyses. 
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The black dots of panel c) in Figure 3 show the correla- 
tion between the electric and magnetic wavelet power as 
a function of kpi. It can be seen that the fluctuations 
are strongly correlated through the inertial range (with 
coefficient ks 1), remain well-correlated between the two 
breakpoints kpt £ (0.45, 2) and begin to loose corre- 
lation quickly above the second breakpoint. A wavelet 
cross-spectral analysis (between 6E y and SB Z was also 
computed; the blue bars show the cross-spectral coher- 
ence, with 1 sigma error bars, also as a function of kpi. 
Again, 6E y and 5B Z are strongly coherent through the 
inertial range and past the first breakpoint. We conclude 
that the electric and magnetic spectra physical and well- 
correlated up to the second spectral breakpoint. Above 
kpi « 2.5 it is difficult to assess the quality of the data. If 
electrostatic waves are present, there is no expectation of 
correlation with SB; however in this initial study, we can- 
not eliminate the possibility of systematic noise at these 
frequencies. Additionally, the effects of low pass filters 
on both the EFW and FGM experiments may modify 
the spectra at these highest (kpi > 3) frequencies. 

To estimate the phase speed of the fluctuations, we use 
Faraday's law and compute the ratio of the electric and 
magnetic spectra. Since the electric field measurements 
are made in the spacecraft (unprimed) frame, we need 
to Lorentz transform to the plasma (primed) frame by 
E = E' + v sw x B. Panel b) of Figure 3 shows the phase 
speed 



v<t,{kpi) 



SE' y {k Pi ) = SEyjkpi) 
SB z (k Pi ) SB z (k Pl ) 



shown as a horizontal bar. Over, and even below, the 
inertial range kpi G (0.015, 0.4) the phase speed is con- 
sistent with the local Alfven speed; this is strong evidence 
of the Alfvcnic nature of the cascade. The red curve in 
panel b) is a fit of the function v a (l + k 2 p 2 ) to the FFT 
curve, where vq is a free parameter which finds a best fit 
at vq ks 55 km/s. This function approximates the disper- 
sion relation of kinetic Alfven waves. The cold-plasma 
whistler wave phase speed goes as v<j, « (kpi) (3~ x / 2 va 
above uj > n c ,-, i.e. linear with kpi, and would form 
a much shallower dispersion above kpi « 1 than is ob- 
served in panel b) of Figure 3. This leads us to believe 
that the Alfven waves in the inertial subrange eventu- 
ally disperse as 'kinetic' Alfven waves above kpi w 1, 
becoming more electrostatic and eventually damping on 
the thermal plasma. Plasma heating by linear dissipa- 
tion of kinetic Alfven waves at/3 ~ 1 has been studied in 
the context of accretion flows [5J. There it was found that 
Landau and transit-time damping contribute to both pro- 
ton and electron heating at short wavelengths, which is 
enhanced for higher (3. Kellogg computed the level of 
electric field fluctuations required to stochastically ther- 
malize protons to 1 AU in the solar wind; he found that 
a spectral density of E 2 w 10~ 11 (V/m) 2 Hz~ 1 was suf- 
ficient. This is one order of magnitude less than our 
observed levels (Figure 2). It is, therefore, plausible to 
conclude that the observed electric spectrum is respon- 
sible for isotropizing the solar wind protons and may be 
the mechanism by which the solar wind maintains its 
fluid-like characteristics. 



(1) 



where v x , v z , B x , and B z are the average x and z compo- 
nents of the solar wind velocity and magnetic field. The 
black dots in panel b) are computed from the wavelet 
spectrum, while the blue line is computed from the FFT 
spectrum. The average Alfven speed va ~ 40 km/s is 
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FIG. 1: Wavelet and time series data of solar wind turbu- 
lence. From the top-down, the five panels show a) the wavelet 
spectrogram of E y , as a function of kpi, b) a similar wavelet 
spectrogram of B z , c) the X and Y components of the mea- 
sured electric field, d) the vector magnetic field, e) plasma 
ion density, f) plasma ion ft, and g) the Alfven Mach number. 
This entire interval was used for the spectral analysis of E y 
and B z . The spectral breakpoints called out. 



FIG. 2: Power spectral density of electric 8E y and magnetic 
fluctuations 5B Z as a function of frequency, computed from 
FFT (black) and Morlet wavelet (red) algorithms. The FFT 
spectrum of electric field (upper panel) shows the effect of 
notch filters and residual spin-tone data. 



FIG. 3: The wavelet (upper) and FFT (lower) power spec- 
tra of E y (green) and B z (black) binned as a function of 
wavenumber kpi (and offset for clarity) in panel a) . The elec- 
tric are multiplied by factor to lie atop the magnetic spec- 
tra. The spectrum is Kolmogorov fc~ 5 / 3 over the interval 
kpi G (0.015,0.45); a spectral breakpoint occurs for both E y 
and B z at kpi w 0.45. A second breakpoint occurs for elec- 
tric spectrum at kpi « 2.5 above which the electric spectrum 
is more exponential. Panel b) shows the ratio of the electric 
to magnetic spectra in the plasma frame; the average Alfven 
speed (va ~ 40 km/s) is shown as a horizontal line. The red 
line is a fitted dispersion curve, discussed in the text. Panel 
c) shows both the cross-coherence of 5E y with 5B Z (as blue 
dots with error bars) and the correlation between the electric 
and magnetic power (as black dots). 
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